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Abstract. We present a short critical overview of different microscopic models for nonrelativistic and 
relativistic magnetoelectric coupling including the so-called "spin current scenario", ab-initio calculations, 
and several recent microscopic approaches to a spin-dependent electric polarization in 3d oxides. 

PACS. 77.80.-e Ferroelectricity and antiferroelectricity - 75.80.+q Magnetomechanical and magnetoelec- 
tric effects, magnetostriction - 71.70.Gm Exchange interactions 
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1 Introduction 

Since Astrov's discovery of the magnetoelectric (ME) ef- 
fect in C^OapQ several microscopic mechanisms of mag- 
netoelectric coupling were proposed]!], however, the mul- 
tiferroicity (see e.g., Refs.[3] and review articles Refs.[21ll]) 
has generated an impressive revival of the activity in the 
field. The microscopic origin of the magnetically driven 
electric polarization is the topic of an intense and con- 
troversial debate. Currently two essentially different spin 
structures of net electric polarization in crystals are con- 
sidered: a bilinear nonrelativistic symmetric spin coupling 



p s = £n mn (s m -s n ) 



(i) 



and a bilinear relativistic antisymmetric spin coupling 

p q = ]TnL[s m x s„] , (2) 

mn 

respectively. If the first term stems somehow or other from 
a spin isotropic Heisenberg exchange interaction, the sec- 
ond term does from antisymmetric Dzyaloshinsky-Moriya 
(DM) coupling. These effective spin-operator forms do 
not discriminate the "ionic" and "electronic" contribu- 
tions to magnetically driven ferroelectricity related with 
the off-center ionic displacements and the electron den- 
sity redistribution, respectively. A microscopic quantum 
theory of the ME effect has not yet been fully developed, 
although several scenarios for particular materials have 
been proposed. Many authors consider that the giant mul- 
tiferroicity requires the existence of sizeable atomic dis- 
placements and structural distortions driven by isotropic 
symmetric exchange coupling [HOSE] or antisymmetric DM 
coupling 8j ignoring, however, the fact that the effects of 



nuclear displacements and electron polarization should be 
described on equal footing, e.g., in frames of the well- 
known shell model of Dick and Overhauser 9 widely used 
in lattice dynamics. Shell and core displacements may be 
of a comparable magnitude. The conventional shell model 
does not take into account the spin and orbital degrees of 
freedom, hence it cannot describe the multiferroic effects. 
In fact, the displacements of both the atomic core and 
electron shell would depend on the spin surroundings pro- 
ducing an synergetic effect of spin-dependent electric po- 
larization. Obviously, this effect manifest itself differently 
in neutron and x-ray diffraction experiments. Sorting out 
two contributions is a key issue in the field. 

The second harmonic generation (SHG) spectroscopy 
reveals a giant effect in compounds with magnetically driven 
ferroelectricity (TbM^Os, MnWCU) thus pointing to an 
electronic rather than ionic origin of the spontaneous polarization|10j 
Interestingly, the exchange-induced electric dipole moment 
([T]), derived many years ago [H](see, also Ref.[12]h seems 
to be a natural electronic mechanism for giant multifer- 
roicity. However, the second, or exchange-relativistic elec- 
tric dipole moment ([2]), despite its visible weakness, is at 
present frequently addressed to be the leading contributor 
to electronic multiferroicity, mainly due to the so-called 
" spin-current" mechanism [131114] . To a large extent it is 
explained by two reasons. First, in a most part of multifer- 
roics the electric polarization is observed when the mag- 
netic ordering is of a type that breaks chiral symmetry, 
e.g., spiral or helical order with a nonzero "spin-current" 
oc [S m x S n ]. Second, such a mechanism allows us to eas- 
ily predict the direction of the ferroelectric polarization for 
certain helical spin structures: P a oc J2mn [Rmn x [Sm x S n ]] 



OC 



mn i 

eijfe(Rmn)fc) [T3], or P a oc [e x Q] (e is the he- 



lix spin rotation axis, Q is the helix wave vector) [14] . It is 
worth noting that phenomenological theory by Mostovoy [Tf 
implies a coexistence of two independent order parameters 
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P(r) and M(r) coupled by the magnetoelectric energy: 
<2> me (P(r), M(r)) = 7 P [M(V • M) - (M • V)M], while for 
the multiferroics under consideration we deal with M(r) 
to be a leading, or first order parameter, and P(r) as a 
second order parameter which is expressed in terms of 
the former one. Furthermore, his phenomenological the- 
ory does not explain neither the origin nor the magnitude 
of the ME coupling parameter 7. 

"Ferroelectricity caused by spin-currents" has estab- 
lished itself as one of the leading paradigms for both the- 
oretical and experimental investigations in the field of 
strong multiferroic coupling. However, a "rule" that chi- 
ral symmetry needs to be broken in order to induce a 
ferroelectric moment at a magnetic phase transition is 
questionable [15 . There are notable exceptions, in partic- 
ular, the manganites RM^Os, HoMn03, where a ferro- 
electric polarization can appear without any indication 
of a magnetic chiral symmetry breaking [5lll6j , and de- 
lafossite CuFei_ x A1^02, where the helimagnetic ordering 
generates a spontaneous electric polarization parallel to 
the helical axis [17] , in sharp contrast with the prediction 
of the spin current model. Moreover, there are increasing 
doubts whether weak exchange-relativistic coupling can 
generate the giant electric polarization observed in some 
multiferroics. 

Another point of hot debates around the microscopic 
origin of ME coupling is related with recent observations of 
a multiferroic behaviour concomitant the incommensurate 
spin spiral ordering in chain cuprates LiVCuC>4 [18,19, 
1213] and LiCu 2 2 [2TU22U25] . At first sight, these cuprates 
seem to be prototypical examples of ID spiral-magnetic 
ferroelectrics revealing the relativistic mechanism of "fer- 
roelectricity caused by spin-currents" [13] • However, these, 
in particular, LiCu2 02 show up a behavior which is ob- 
viously counterintuitive within the framework of spiral- 
magnetic ferroelectricity [5T]. Furthermore, quantum he- 
limagnets NaCu202 and Li 2 ZrCu04 with a very similar 
C11O2 spin chain arrangement do not reveal signatures 
of a multiferroic behavior. Thus, there is no clear under- 
standing connecting all these striking properties. 

Despite its popularity, the original " spin-current" model 
and its later versions [2"4"U2"5] (see also Ref. [SB]) seem to 
be questionable as the authors proceed with an unrealis- 
tic scenario. Indeed, when addressing a generic centrosym- 
metric M1-O-M2 system they groundlessly assume an ef- 
fective spin polarization Zeeman field ( m i ^ s a l° ca l 
magnetic moment) to align noncollinearly the spins of 3d 
electrons and to provide a nonzero value of the two-site 
spin current [Si x 82]- The energy separation U originates 
from the local Coulomb repulsion and the Hund coupling 
in the magnetically ordered phase 13,24,25,26 . Such an 
assumption goes beyond all the thinkable perturbation 
schemes and leads to an unphysically large effect of break- 
ing of a spatial symmetry induced by a spin configuration 
that manifests itself in an emergence of a nonzero electric 
dipole moment for an isolated centrosymmetric M1-O-M2 
system [13] . 

Size of the macroscopic polarization P in nonmagnetic 
ferroelectrics computed by modern ab-initio band struc- 



ture methods agrees exceptionally well with the ones ob- 
served experimentally. However, the state of the art ab- 
initio computations for different multiferroics: mangan- 
ites HoMn0 3 [2I], TbMn 2 5 [28], HoMn 2 5 [29], spin spi- 
ral chain cuprates L1VCUO4 and LiCu202 [30] yield data 
spread within one-two orders of magnitude with ambigu- 
ous and unreasonable values of polarization depending 
on whether these make use of theoretical or experimen- 
tal structural data or different values of the correlation 
parameters. The basic starting points of the current ver- 
sions of such spin-polarized approaches as the LSDA seem 
to exclude any possibility to obtain a reliable quantitative 
estimation of the spin-dependent electric polarization in 
multiferroics [31] . We should emphasize two weak points 
of so-called first-principle calculations which appear as 
usual to be well forgotten in the literature. First, these 
approaches imply the spin configuration induces immedi- 
ately the appropriate breaking of spatial symmetry that 
makes the symmetry-breaking effect of a spin configu- 
ration to be unphysically large. It is worth noting that 
the spin-current scenario [13] starts with the same LSDA- 
like assumption of unphysically large symmetry-breaking 
spin-magnetic field. Conventional schemes imply just the 
opposite, however, a physically more reasonable picture 
when the charge and orbital anisotropics induce a spin 
anisotropy. Second, the first-principle calculations neglect 
quantum fluctuations, that restricts drastically their ap- 
plicability to a correct description of the ME coupling de- 
rived from the high-order perturbation effects. 

Below we present a short overview of different micro- 
scopic approaches to spin-dependent electric polarization. 
In Sec. II we address a systematic standard microscopic 
theory which implies the derivation of effective spin op- 
erators for nonrelativistic and relativistic contributions to 
electric polarization of the generic three-site two-hole clus- 
ter such as Cui-0-Cu2. In Sec. Ill we address an alter- 
native approach based on the parity breaking exchange 
coupling and the exchange-induced electric polarization 
effects. Sec. IV is focused on the microscopic origin of the 
multiferroic behaviour observed in the edge-shared Cu02 
chain compounds LiVCu04 and LiCu202- In Sec.V we 
Afyaw attention to anomalous magnetoelectric properties 
of the electron-hole dimers to be precursors of the dis- 
proportionated phase which droplets can survive even in 
nominally pure undoped manganites. 



2 Spin-dependent electric polarization 
three-site Mi-0-M 2 cluster 



in a 



Generic three-site M1-O-M2 cluster forms a basic element 
of the crystalline and electron structure of 3d oxides. A 
realistic perturbation scheme needed to describe properly 
the active M 3d and O 2p electron states implies strong 
intra-atomic correlations, the comparable effect of crystal 
field, the quenching of orbital moments by a low-symmetry 
crystal field, account for the dp-transfer up to the fourth 
order effects, and a rather small spin-orbital coupling. To 
this end we make use of a technique suggested in refs. [351 
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2.1 Three-site two-hole Mi-0-M 2 cluster 

For illustration, below we address a three-site (Cui-O- 
CU2) two-hole system typical for cuprates with a tetrago- 
nal Cu on-site symmetry and a Cu 3d x 2_ y 2 ground states. 
We start with the construction of the spin-singlet and 
spin-triplet wave functions for the system taking account 
of the p-d hopping, on-site hole-hole repulsion, and crystal 
field effects for excited configurations {n} (Oil, 110, 020, 
200, 002) with different hole occupation of Cui, O, and 
CU2 sites, respectively. The p-d hopping for Cu-0 bond 
implies a conventional Hamiltonian 

Hpd = Y tpetdpPadp + h.C. , (3) 

a/3 

where creates a hole in the a state on the oxygen site, 
while dp annihilates a hole in the (3 state on the copper 
site; tpadp is a pd-transfer integral. For basic 101 configu- 
ration with two d x 2_ y 2 holes localized on its parent sites 
we arrive at a perturbed wave function as follows 

&W\;SM = Vs[$l0l;SM + ^"it^WJSAf], 

r{n}^101 

(4) 

where the summation runs both on different configura- 
tions and different orbital r states; 

vs = [i+J2 \c {n} (sr)\T 1/2 (5) 

{n}r 

is a normalization factor. It is worth noting that the prob- 
ability amplitudes, or hybridization parameters, c/qh}, 
C{no} °c tpd, C{200}, C{ 20}, C{ 02} oc t 2 pd . To account for 
relativistic effects in the three-site cluster one should in- 
corporate the spin-orbital coupling V so = Yli£nl{h ■ &i) 
both for 3d- and 2p-holes with a single particle constant 

> for electrons and < for holes. 

In terms of the hole spins the conventional bilinear spin 
Hamiltonian for the CU1-O-CU2 system reads as follows: 

#.(12) = J(si • s 2 ) + D ■ [si x s 2 ] + S!Ks 2 , (6) 

where J oc t^ d is an exchange integral, D oc t pd ^ n i is a 

Dzyaloshinsky vector, K oc ip d £ 2 ; is a symmetric second- 
rank tensor of the anisotropy constants [3"2"]l33] . 

2.2 Nonrelativistic mechanism of spin-dependent 
electric polarization:local and nonlocal terms 

Projecting the electric dipole moment P = |e|(ri + r 2 ) 
on the spin singlet or triplet ground state of the two-hole 
system we arrive at an effective electric polarization of 
the three-center system (P) s = (^ioi ; sm|P|*ioi ; sm) to 
consist of local and nonlocal terms, which accomodate 
the diagonal and off-diagonal on the ionic configurations 
matrix elements, respectively 31 . The local contribution 



3 

describes the redistribution of the local on-site charge den- 
sity and can be written as follows: 

(P) l s ° cal = \e\\ Vs \ 2 [(Ri + R 2 + (Ri + Ro) £ |cno(Sr)| 2 

r 

+ (Ro + R 2 I con^r) 1 2 + 2R G |c 02 o(Sr)| 2 
r r 

+ 2R 1 £| C20 o(^)| 2 + 2R2£|coo2(5r)| 2 ] -P , (7) 
r r 

where Po = |e|(Ri + R2) is a bare purely ionic two- hole 
dipole moment. Obviously, the net local electric polariza- 
tion depends only on R.y vectors (Rio, R20, Ri2)- It is 
worth noting that the net local electric polarization lies 
in the Cui-0-Cu2 plane. The nonlocal, or overlap con- 
tribution is related with the off-diagonal two-site matrix 
elements of P [31] . 

The effective electric polarization differs for the singlet 
and triplet pairing due to a respective singlet-triplet differ- 
ence in the hybridization amplitudes C{ n y(Sr). Hence we 
may introduce an effective nonrelativistic exchange- dipole 
spin operator 

p s = n(§! ■ § 2 ) (8) 

with an exchange- dipole moment 

n = <p) t -(p) s , (9) 

which can be easily deduced from Exp. ([7]). The effec- 
tive nonrelativistic exchange-dipolc moment is determined 
by competitive local and nonlocal contributions of sev- 
eral configurations [31] . It is worth noting that for the 
collinear CU1-O-CU2 bonding both contributions vanish. 
As a whole, the exchange-dipole moment vanishes, if the 
M1-O-M2 cluster has a center of symmetry. 

It is worth noting that we addressed only the charge 
density redistribution effects for Cu 3d and O 2p states 
and neglected a direct electronic polarization effects for 
the both metal and anion ions. These effects be incorpo- 
rated to the theory, if other orbitals, e.g. ns- for oxygen 
ion, will be included in the initial orbital basis set. To 
proceed with these effects an alternative approach may be 
applied by using a generalized shell model [34], 



2.3 Relativistic mechanism of the spin-dependent 
electric polarization 

At variance with a scenario by Katsura et al. |13j we have 
applied a conventional procedure to derive an effective 
spin operator for a relativistic contribution to the elec- 
tric dipole moment in the three-site M1-O-M2 system like 
a technique suggested in references [32 , 33 to derive ex- 
pressions for the Cu and O spin-orbital contributions to 
the DM coupling in cuprates. 

The spin-orbital coupling Vso for copper and oxygen 
ions drives the singlet-triplet mixing which gives rise to 
a relativistic contribution to electric polarization deduced 
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from an effective spin operator, or an exchange-relativistic- 
dipole moment 



-irr = n[si x § 2 ] , 



(10) 



where 



n, 



-i{* s \PiWti) = m\m$s 



(# t |Pii#t»^ 



is an exchange-relativistic-dipole tensor (iP s and W t j are 
spin singlet and spin triplet wave functions ((3]), respec- 
tively). It is easy to see that this quantity has a clear 
physical meaning to be in fact a dipole matrix element for 
a singlet-triplet electro-dipole transition in our three-site 
cluster. Taking into account equation ©, we arrive at a 
simple form for the exchange-relativistic-dipole moment 
as follows 

P = -in(D-[si xs 2 ]) . (12) 

It is worth noting that this vector lies in the Cui-0-Cu2 
plane and its direction, at variance with the spin-current 
model [13] , does not depend on the spin configuration. Fur- 
thermore, we see that the both nonrelativistic and rel- 
ativistic contributions to effective dipole moment in the 
Cui-0-Cu2 system have the same direction: P Sja oc nf 2 - 
In other words, in the both cases the spin-correlation fac- 
tors, (§i • §2) and [si x §2], do modulate a pre-existing 
dipole moment. The DM type exchange-relativistic-dipole 
moment (|12j) is believed to be a dominant relativistic con- 
tribution to the electric polarization in a CU1-O-CU2 clus- 
ter. It is worth noting that the exchange-dipole moment 
operator ([5]) and exchange-relativistic-dipole moment op- 
erator (fl~2"|) are obvious counterparts of the Heisenberg 
symmetric exchange and DM antisymmetric exchange, re- 
spectively. Hence, the Moriya like relation [55] |77y| ~ 
Ag/g\U\ seems to be a reasonable estimation for the re- 
sultant relativistic contribution to electric polarization in 
M1-O-M2 clusters. At present, it is a difficult and, proba- 
bly, hopeless task to propose a more reliable and so physi- 
cally clear estimate. Taking a typical value of Ag/g ~ 0.1 
we estimate the maximal value of |iTy| as 10 _3 |e|A(~ 
10 2 [iC I to 2 ) that points to the exchange-relativistic mech- 
anism to be a weak contributor to a giant multiferroicity 
with ferroelectric polarization of the order of 10 3 /iC/m 2 as 
in TbMnC>3 [3] , though it may be a noticeable contributor 
in, e.g., Ni 3 V 2 8 [Mi- 



3 Parity breaking exchange coupling and 
exchange-induced electric polarization 

Along with many advantages of the three-site cluster model 
it has a clear imperfection not uncovering a direct role 
played by exchange coupling as a driving force to induce 
a spin-dependent electric polarization. Below we'll address 
an alternative approach starting with a spin center such 
as a MeO n cluster in 3d oxides exchange-coupled with 
a magnetic surroundings. Then the magnetoelectric cou- 
pling can be related to the spin-dependent electric fields 



generated by the spin surroundings in a magnetic crystal. 
In this connection we should point out some properties 
of exchange interaction that usually are missed in conven- 
tional treatment of Heisenberg exchange coupling. Follow- 
ing the paper by Tanabe et a/.[TT] (see, also Ref . [T2"] ) we 
do start with a simple introduction to exchange-induced 
electric polarization effects. 

Let address the one-particle (electron/hole) center in 
a crystallographically centrosymmetric position of a mag- 
netic crystal. Then all the particle states can be of a cer- 
tain spatial parity, even (g) or odd (u), respectively. Hav- 
ing in mind the 3d centers we'll assume an even-parity 
ground state \g). For simplicity we restrict ourselves to 
only one excited odd-parity state \u). The exchange cou- 
pling with the surrounding spins can be written as follows: 



V,: 



= £)j(R n )(s-S n ), 



where I(R n ) is an orbital operator with a matrix 



^99 (^"m) ^gu (^n) 
^-ug (-^n) J-uu (f^n) 



(13) 



(14) 



The parity-breaking off-diagonal part of exchange cou- 
pling can lift the center of symmetry and mix \g) and 
\u) states thus resulting in a nonzero electric dipole polar- 
ization of the ground state 



, H n (s • S ra ) . 



where 



with A 



Fin — 27g U (R, n ) 



(ff|er| 



A 



(15) 



(16) 



11 9 



I'!-) 



u — t g . It is easy to see that in frames of 
a mean-field approximation the nonzero dipole moment 
shows up only for spin-noncentrosymmetric surroundings, 
that is if the condition (S(R„)) = (S(— R„)) is broken. 
For an isotropic bilinear exchange coupling this implies a 
spin frustration. 

It should be noted that at variance with the spin- 
current model [13] the direction of the exchange-induced 
dipole moment for i,j pair does not depend on the direc- 
tion of spins Si and Sj. In other words, the spin-correlation 
factor (Si ■ Sj) modulates a pre-existing dipole moment Tl 
which direction and value depend on the Mei-O-Mej bond 
geometry and orbitals involved in the exchange coupling. 

The magnitude of the off-diagonal exchange integrals 
can sufficiently exceed that of a conventional diagonal ex- 
change integral. Given reasonable estimations for the off- 
diagonal exchange integrals I ug « 0.1 eV, the u — g en- 
ergy separation A ug « 2 eV, the dipole matrix element 
|(g|er|u)| « O.lA, spin function |((s • S n ))| w 1 we ar- 
rive at an estimation of the maximal value of the elec- 
tric polarization: P ps 10 4 /iC/m 2 . This estimate points 
to the exchange-induced electric polarization to be po- 
tentially the most significant source of magnetoelectric 
coupling for new giant multiferroics. It is worth noting 
that the exchange-induced polarization effect we consider 
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is particularly strong for the 3d clusters such as MeO„ 
with an intensive low-lying electro-dipole allowed tran- 
sition \g) — > \u) which initial and final states are cou- 
pled due to a strong exchange interaction with a spin sur- 
roundings [37] . This simple rule may be practically used 
to search for new multiferroic materials. 

The parity-breaking exchange coupling can produce a 
strong electric polarization of oxygen ions in 3d oxides 
which can be written as follows 

p = ^n n ((s )-s„), (17) 

n 

where S n are the spins of the surrounding 3d ions, (So) oc 

Y] I n S„ is a spin polarization of the oxygen ion due to 

the surrounding 3d ions with I „ being an exchange cou- 
pling tensor. It seems the oxygen exchange-induced elec- 
tric polarization of purely electron origin has been too 
little appreciated in the current pictures of multiferroicity 
in 3d oxides. 



4 Origin of multiferroic properties in the 
edge-shared Cu0 2 chain compounds 

4.1 Cancellation of the spin-dependent electric 
polarization in perfect edge-shared CuOj chains 

Recent observations of a multiferroic behaviour concomi- 
tant the incommensurate spin spiral ordering in cuprates 
LiVCu0 4 p l ITI p Ol and LiCu 2 2 [21] with nearly per- 
fect edge-shared CuC>2 chains (see Fig.rjJ) challenge the 
multiferroic community. From the viewpoint of the spin- 
current model, these cuprates seem to be prototypical 
examples of the ID spiral-magnetic ferroelectrics reveal- 
ing the relativistic mechanism of " ferroelectricity caused 
by spin-currents" [13]. However, as we see from discus- 
sion above, isolated perfect edge-shared centrosymmetric 
CuC>2 chains cannot produce a spin-dependent electric po- 
larization both of nonrelativistic and relativistic origin. In- 
deed, the net nonrelativistic polarization of a spin chain 
formed by Cu ions positioned at the center of symmetry 
can be written as follows [11] 

P e// = n ]T [{Sj ■ S j+1 ) - (3 r S,-_r)] , (18) 

j=even 

hence for a simple plane spiral ordering in perfect edge- 
shared CuC>2 chains we arrive at a twofold cancellation 
effect due to the zeroth value both of the II and the spin- 
correllation factor in brackets. A twofold cancellation ef- 
fect takes place for the relativistic contribution (fl~2")) to the 
spin-dependent electric polarization as well, because both 
the exchange-dipole moment II and Dzyaloshinsky vector 
D turn into zero. Indeed, a specific symmetry of Cui-O- 
Cu2 bonds in edge-shared CuC>2 chains (see Fig.[T]) results 
in a full cancellation of a net Dzyaloshinsky vector, though 
the partial Cui-0/ j //-Cu2 contributions survive being of 
opposite signs 32,33"]. 




Fig. 1. (Color online) The fragment of a typical edge- 
shared CuC>2 chain. Note the antiparallel orientation of the 
Dzyaloshinsky vectors in the Cui-Oj-Cu2 and Cui-Oij-Cu2 
bonds both directed perpendicular to the chain plane 32 . 

The absence of a spin-dependent ferroelectric polariza- 
tion in perfect edge-shared CuC>2 chains is a simple corol- 
lary of its centrosymmetry. The nonzero effect predicted 
by the spin-current model [13] is related to an unphysical 
symmetry breaking engendered by a strong Active "exter- 
nal" nonuniform field needed to align spirally the chain 
spins. Thus we may state that the edge-shared CuC>4 pla- 
qucttes arrangement in the CuC>2 chains appears to be 
robust regarding the inducing of the spin-dependent elec- 
tric polarization both of the nonrelativistic and relativistic 
origin. It means that we should look for the origin of puz- 
zling multiferroicity observed in L1VCUO4 and LiCu202 
somewhere within the out-of-chain stuff. 



4.2 Nonstoichiometry and multiferroic behaviour of 
edge-sharing CuOo chain compounds LiVCuO j and 
LiCu 2 2 

According to the spin-current theory [13] a net electric po- 
larization induced by a spin-spiral ordering in CuC>2 chains 
of quantum helimagnets LiVCuO^ , NaCu2C>2, LiCu2C>2 , 
and Li2ZrCuC>4 with a very similar CuC>2 spin chain ar- 
rangement is directed as shown in Fig.[T]with a magnitude 
proportional to sin where <fi is the pitch angle. Thus we 
should anticipate comparable values of a net chain electric 
polarization in these cuprates with pitch angles 85°, 82°, 
62°, 33°, respectively 38J. However, these cuprates show 
up a behavior which cannot be explained within the frame- 
work of spiral-magnetic ferroelectricity [131114] . First, in 
accordance with a cancellation rule discussed above the 
quantum helimagnets NaCu20 2 and Li 2 ZrCu04 do not 
reveal any signatures of a multiferroic behavior while the 
both LiVCu04 and LiCu2C>2 systems reveal a mysteri- 
ous behavior with conflicting results obtained by different 
groups. Indeed, Yasui et al.^9\ claim that LiVCuC>4 re- 
veals clear deviations from the predictions of spin-current 
models [14, 13] while Schrettle et aL[20] assure of its ap- 
plicability. In contrast to LiVCuC>4, LICU2O2 shows up 
a behavior which is obviously counterintuitive within the 
framework of spiral-magnetic ferroelectricity [21] (see Fig. [2]) 
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Fig. 2. (Color online) Direction of ferroelectric polarization in 
LiCu2 02 for different spin spiral plane orientation as observed 
by Park et al. [21] and predicted by the " nonstoichiometry" 
mechanism [40] . For comparison the predictions of the spin- 
current model [13J.14, are shown. 



It is worth noting that at variance with Park et aL[21], 
Naito et al. [18] have not found any evidence for ferroelec- 
tric anomalies in L1CU2O2. The ferroelectric anomaly in 
LiVCu04 reveals a magnitude (P a ~ 30/iC/m 2 ) compara- 
ble to that of the multiferroic N13V2O8 [36 while LiCu202 
shows up an order of magnitude lesser effect [21] . Recently 
we have shown that the unconventional multiferroic be- 
haviour observed in samples of L1VCUO4 and LiCu2 02 
can have nothing to do with " spin-currents" , this can be 
related with a nonstoichiometry in samples under consid- 
eration [39,40 . Their " multiferroicity" can be consistently 
explained if one takes into account the nonrelativistic 
exchange-induced electric polarization on the Cu 2+ cen- 
ters substituting for the positions native for Li-ions in 
LiVCu04 and Cu 1+ -ions in LiCu202, respectively [3"51l4"0] . 
Such a mechanism does explain even subtle features of a 
multiferroic behavior observed in L1VCUO4 and LiCu 2 02. 
These results raise a number of questions of great impor- 
tance for physics of magnetism and multiferroicity in the 
spin s=l/2 quantum matter of LiVCuCU and LiCu20 2 . 
Whether a multiferroic behavior would be observed in sto- 
ichiometric samples with a regular arrangement of Cu 2+ 
and Cu 1+ ions? Recently, in order to exclude the non- 
stoichiometry as a source of a multiferroic behavior, the 
Nagoya University group has prepared single-crystal sam- 
ples of LiCu202 with a controlled stoichiometry, which 
had "neither the atomic deficiency nor the mixing of Cu 
and Li atoms" . The authors have tried to detect an electric 
polarization, however, at variance with earlier findings [21i 
22 , 25j have not found any systematic data [U] . Hardly vis- 
ible bumps of the capacitance for these samples {AC /C « 
0.001 at E || c) observed at critical temperatures T^ and 
Tjv 2 @2] cannot validate the spin-current theory [T^UTl] . 
These results are believed to support strongly the cancel- 
lation rule at work in the edge-shared Cu02 chains and the 
nonstoichiometry as a source of a multiferroic behavior ob- 
served earlier both in LiVCu0 4 PMD] and LiCu 2 2 [22 

1220231 • 



5 Electron-hole dimers, electronic phase 
separation, and dielectric anomalies in 
undoped parent manganites 

Even the nominally pure globally centrosymmetric par- 
ent manganite LaMn03 exhibits a puzzling multiferroic- 
like behavior inconsistent with a simple picture of an A- 
type antifcrromagnetic insulator (A-AFI) with a cooper- 
ative Jahn- Teller ordering. Its anomalous properties are 
assigned to charge transfer instabilities and competition 
between insulating A-AFI phase and metallic-like dynam- 
ically disproportionated phase formally separated by a 
first-order phase transition at Tdisp = T jt^ 750 K [4"3"] . 
The unconventional high-temperature phase is addressed 
to be a specific electron-hole Bose liquid (EHBL) rather 
than a simple "chemically" disproportionated R(Mn 2+ Mn 4+ )03 
phase. New phase does nucleate as a result of the charge 
transfer (CT) instability and evolves from the self-trapped 
CT excitons, or specific EH-dimers, which seem to be a 
precursor of both insulating and metallic-like ferromag- 
netic phases observed in manganites. The view of a self- 
trapped CT exciton to model a Mn 2+ -Mn 4+ pair is typical 
for a chemical view of disproportionation, and is strongly 
oversimplified. Actually we deal with an EH-dimer to be 
a dynamically charge fluctuating system of coupled elec- 
tron MnOg 0- and hole MnO^ - centers having been glued 
in a lattice due to strong electron-lattice polarization ef- 
fects. In other words, we should proceed with a rather 
complex physical view of disproportionation phenomena 
which first implies a charge exchange reaction 

Mn 2+ + Mn 4+ <-► Mn 4+ + Mn 2+ , (19) 

governed by a two-particle charge transfer integral 

t B = (Mn 2+ Mn 4+ | J tf s |Mn 4+ Mn 2+ ) , (20) 

where Hb is an effective two-particle (bosonic) transfer 
Hamiltonian, and we assume a parallel orientation of all 
the spins. As a result of this quantum process the bare 
ionic states with site-centred charge order and the same 
bare energy Eq transform into two EH-dimer states with 
an indefinite valence and bond-centred charge-order 

|±) = -^(|Mn 2+ Mn 4+ ) ± |Mn 4+ Mn 2+ )) (21) 
v2 

with the energies E± = Eq ± ts- In other words, the 
exchange reaction restores the bare charge symmetry. In 
both |±) states the site manganese valence is indefinite 
with quantum fluctuations between +2 and +4, however, 
with a mean value of +3. Interestingly, in contrast with the 
ionic states, the EH-dimer states |±) have both a distinct 
electron-hole and inversion symmetry, even parity (s-type 
symmetry) for |+), and odd parity (p-type symmetry) for 
|— ) states, respectively. Both states are coupled by a large 
electric-dipole matrix element: 

(+|d|-) = 2eR Mn Mn, (22) 
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Fig. 3. (Color online) Spin structure of the EH-dimer with a 
step-by-step inclusion of one- and two-particle charge transfer. 
Arrows point to electric dipole moment for bare site-centred 
dimer configurations. 



First of all, one should note the relatively high static di- 
electric constant in LaMnC>3 at T = (eq ~ 18 — 20) ap- 
proaching to values typical for genuine multifcrroic sys- 
tems (t?o ~ 25), whereas for the conventional nonpolar sys- 
tems, Eq varies within 1-5. The entire e'(lo,T)- T pattern 
across 77-900 T has two prominent features: (i) near TV 
and (ii) near Tjt to be essential signatures of puzzlingly 
unexpected multiferroicity, however, the intrinsic electri- 
cal polarization probably develops locally with no global 
ferroelectric order. The observation of an intrinsic dielec- 
tric response in the globally centrosymmetric LaMn03, 
where no ferroelectric order is possible due to the ab- 
sence of off-centre distortion in MnOe octahedra cannot 
be explained in frames of a conventional uniform anti- 
ferromagnetic insulating A-AFI scenario and agrees with 
the electronic A-AFI/EHBL phase separated state with 
a coexistence of the non-polar A-AFI phase and a highly 
polarizablc EHBL phase |43J. 



where RmuMu is the Mn-Mn separation. In a nonrelativis- 
tic approximation the spin structure of the EH-dimer will 
be determined by the isotropic Heisenberg exchange cou- 
pling V ex = J (Si • S 2 ), and the two-particle charge trans- 
fer characterized by a respective transfer integral which 
depend on the spin states. Both terms can be easily diag- 
onalized in the net spin S representation so that for the 
energy we arrive at 

/ 25 1 

E S = -[S(S + 1) - -] ± -S(S + 1) t B , (23) 

where ± corresponds to two quantum superpositions |±) 
with s- and p-type symmetry, respectively. It is worth not- 
ing that the bosonic double exchange contribution for- 
mally corresponds to a ferromagnetic exchange coupling 
with Jb = — Tol^sl- We see that the cumulative effect of 
the Heisenberg exchange and the bosonic double exchange 
results in a stabilization of the S = 4 high-spin (ferromag- 
netic) state of the EH-dimer provided \ts\ > 10 J (see 
Fig[3]) and the S = 1 low-spin (ferrimagnetic) state other- 
wise. The spin states with intermediate S values: S= 2, 3 
correspond to a classical noncollinear ordering. 

The EH-dimer reveals unconventional magnetoelectric 
properties. Indeed, the two-particle bosonic transport and 
respective kinetic contribution to stabilization of the ferro- 
magnetic ordering can be suppressed by a relatively small 
electric fields that makes the EH-dimer to be a promis- 
ing magnetoelectric cell especially for the heavy rare-earth 
manganites RMn03 (R=Dy, Ho, Y, Er) with supposedly a 
ferro-antiferro instability. In addition, a strong anisotropy 
of the dimer's electric polarizability is noteworthy. In an 
external electric field the EH-dimers tend to align along 
the field. 

Anomalous electric polarisability of the EH dimers and 
EH droplets that would result in dielectric anomalies in 
the EHBL phase and the phase-separated state of LaMnC>3 . 
Indeed, such anomalies were reported recently both for 
poly- and single-crystalline samples of the parent LaMn03, z 



6 Conclusion 

We have considered several mechanisms of spin-dependent 
electric polarization in 3d oxides. Starting with a generic 
three-site two-hole cluster and a realistic perturbation scheme 
we have deduced both nonrelativistic and relativistic con- 
tributions to the electric polarization. Nonrelativistic mech- 
anism related to the redistribution of the local on-site 
charge density due to the pd covalency and the exchange 
coupling is believed to govern the multifcrroic behaviour 
in 3d oxides. The approach realized has much in common 
with the mechanism of the bond- and site-centered charge 
order competition (see, e.g. Ref.[15]) though we started 
with the elementary pd charge transfer rather than the 
dd charge transfer. An alternative approach to the deriva- 
tion of the spin-dependent electric polarization was con- 
sidered which is based on the parity-breaking exchange 
coupling and the exchange induced polarization. As an 
actual application of the microscopic approach we discuss 
recent observations of multifcrroic behaviour concomitant 
the incommensurate spin spiral ordering in s=l/2 chain 
cuprates LiVCuC"4 and LiCu202- We argued that the mul- 
tiferroicity observed in these nonstoichiometric cuprate 
samples has nothing to do with " spin currents" and can be 
consistently explained if one takes into account the non- 
relativistic exchange-induced electric polarization on the 
Cu 2+ centers substituting for the positions native for the 
Cu + -ions in LiCu202 or the positions native for the Li + - 
ions in L1VCUO4, respectively. We argued that a charge 
transfer instability accompanied by nucleation of the electron- 
hole dimers and droplets in 3d oxides gives rise to a novel 
type of magnetoelectric coupling due to a field-induced re- 
distribution of the electron-hole droplet volume fraction. 
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